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Abstract. We present detailed three-dimensional mod- 
elling of the stellar luminosity distribution for the disks 
of 31 relatively nearby (< 110 Mpc) edge-on spiral galax- 
ies. In contrast to most of the standard methods avail- 
able in the literature we take into account the full three- 
dimensional information of the disk. We minimize the 
difference between the observed 2D-image and an image 
of our 3D-disk model integrated along the line of sight. 
Thereby we specify the inclination, the fitting function 
for the z-distribution of the disk, and the best values for 
the structural parameters such as scalelength, scaleheight, 
central surface brightness, and a disk cut-off radius. From 
a comparison of two independently developed methods we 
conclude, that the discrepancies e.g. for the scaleheights 
and scalelengths are of the order of « 10%. These dif- 
ferences are not due to the individual method itself, but 
rather to the selected fitting region, which masks the bulge 
component, the dust lane, or present foreground stars. 
Other serious limitations are small but appreciable intrin- 
sic deviations of real disks compared to the simple input 
model. 

In this paper we describe the methods and present contour 
plots as well as radial profiles for all galaxies without pre- 
viously published surface photometry. Resulting parame- 
ters are given for the complete sample. 

Key words: galaxies: fundamental parameters - galaxies: 
surface photometry - galaxies: structure - galaxies: spiral 



1. Introduction 

Global parameters of galactic disks can be used to con- 
strain the formation process as well as the evolution of 
disk galaxies. Recently it has become possible to deduce 
parameters such as scalelength or central surface density 
from numerical or semi-analytical self-consistent galaxy 
formation models (Syer et al. 



1999, Mo et al. 1998, Dal 



canton et al. 1997) and compare the results with observed 



* Based on observations collected at the European South- 
ern Observatory, Chile and Lowell Observatory, Flagstaff (AZ), 
USA 



values. Observationally the Hubble Deep Field gives the 
opportunity to study morphological features and simple 
structural parameters for galaxies even at high redshifts 



1998) 



(Takamiya [1999| , Marleau & Simard |1998| , Fasano & Fil- 
ippi 



Former statistical studies providing sets of homogeneously 
derived parameters for nearby galaxy samples are those of 



de Jong (1996b) using a two-dimensional and of Courteau 
( 1996| ) with a one-dimensional decomposition technique. 
However, so far only a few statistical studies based on 



high qua lity C CD data (de Grijs 1998, Barteldrees & 



Dettmar [1994| , hereafter Paper I) have addressed the ac- 
tual three-dimensional structure of disk galaxies with re- 
gard to the stellar distribution perpendicular to the plane 
(z-direction) and taking into account an outer truncation 
(cut-off radius) for th e disk, first introduced by van der 
Kruit & Searle ( |l981a| ). 

While the stellar distribution perpendicular to the plane 
(z-profile) could result from various "heating" processes 

Hern- 

quist 



during the galactic evolution (Toth 
Mihos 



1995), the cut-off radius can be used to 



constrain either the angular momentum of a protocloud 
(van der Kruit 1987) or possible starformation thresholds 
in the gaseous disk (Kennicutt 1989). 
In the following, we have compiled parameters of galactic 
disks for a total of 31 edge-on galaxies in different optical 
filters of which 17 objects have been already discussed in 
Paper I. One goal of our detailed comparison of several 
independent fitting procedures is to study the influence 
of the applied techniques on resulting disk parameters. It 
also provides the data base for statistical analysis address- 
ing some of the above mentioned questions (e.g., Pohlen 
et al. 200C). For 14 objects surface photometric data are 



given for the first time in Appendix A. 

2. Observations and data reductions 

2.1. Observations 

The observations were carried out at the 42-inch (1.06m) 
telescope of the Lowell Observatory located on the An- 
derson Mesa dark side during several nights in December 
1988 (run identification: LI) and January 1989 (L2) and 
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at the 2.2m telescope at ESO/La Silla during three runs in 
June 1985 (El), March (E2) and June (E3) 1987. At the 
42-inch telescope we used a 2:1 focal reducer with the f/8 
secondary, equipped with a CCD camera which is based 
on a thinned TI 800x800 WFPC 1 CCD with 15 fiui pixel- 
size resulting in a field of approximately 9' with a scale of 
0.7"pixer^ Ima ges were taken with a standard Johnson 
R filter. Observations at the 2.2m telescope were carried 
out with the ESQ CCD adapter using a 512x320 RCA 
chip, giving an effective field size of « 3' x 2' and a scale 
of 0.36"pixel~^. For the ESO observing ru ns we used the 
g, r, and i filters of the Thuan and Gunn ( 1976 ) system. 
Exposures were mainly taken in the g or r band, and only 
seven galaxies were observed in all three filters. 

2.2. Sample selection 

The northern sample observed at Lowell was selected au- 
tomatically in an electronic version of the UGC-catalog 



(Nilson 1973) searching for galaxies with an inclination 
class 7 matching the field size. After visual inspection 
to check the inclination and remove interacting and dis- 
turbed galaxies, the observed sample was chosen out of 
the remaining galaxies according to the allocated observ- 
ing time. 

For the southern sky there is no comparable catalog pro- 
viding information of inclinations directly. Using the axial 
ratios given e.g. in the ESO-Lauberts & Valentijn catalog 
( |1989D will introduce a selection bias preferring late type 
galaxies with lower B/D ratio (Guthrie 1992, Bottinelli et 
al. |I98§). One way to avoid this is extending the first se- 
lection to much lower axis ratios, comparable to i « 65°, 
and then checking the inclination by eye. Therefore we se- 
lected the galaxies according to the field size of about 2' 
from a visual inspection of film copies of the southern sky 
survey (see Paper I). 

Table |] gives a list of the resulting sample used during 
our fitting process, with a serial number (i), the principal 
galaxy name (2)^ the used filter (3)^ the integration time in 
minutes U), and the run identification label (5), whereas 
the marks images already published in Paper I. Fur- 
ther parameters are taken from the RC3 catalogue (de 
Vaucouleurs et al. |1991| ): the right ascension (6) and dec- 
lination ('?;, the RC3 coded Hubble-type fsj, the Hubble 
parameter T (9), and the D25 diameter in arcminutes (10). 
In the case of ESO 578-025 parameters are taken from the 



ESO-Uppsala catalogue (Lauberts 1982) 



2.3. Data reduction 

We applied standard reduction techniques for bias sub- 
traction, bad pixel correction and flatfielding. Following 
the procedure described in Paper I the sky background 
was fitted for each image using the edges of the individ- 
ual frames to reduce any large scale inhomogeneity in the 
field of view. For part of the data we tried to remove the 



foreground stars from the image, but even with sophisti- 
cated PSF fitting using IRAF-DAOphot routines we were 
not able to remove stars without any confusion. The re- 
maining residuals were always of the order of the discussed 
signal. This technique could only be used to mask out the 
regions affected by stars. In order to increase the signal- 
to-noise ratio near the level of the sky background part 
of the data was filtered using a weighted smoothing al- 
gorithm (see Paper I). Thereby the noise was reduced by 
about one magnitude measured with a three sigma devia- 
tion on the background, whereas the interpretation of the 
faint structure is hampered by this process. 
We therefore conclude that the best way is to omit any 
additional image modifications, besides a rotation of the 
image to the major axis of the disk. 



2.4. Photometric calibration 

Most of the images were taken during non-photometric 
nights, therefore we tried a different way to perform pho- 
tometric calibration. Comparing simulated aperture mea- 
surements with published integrated aperture data led to 
the best possible homogeneous calibration of the whole 
sample. Most of the southern galaxies wer e calib rated us- 
ing the catalogue of Lauberts & Valentijn (1989), whereas 
nedQ was used for all northern galaxies. We used equa- 
tion (1) derived in Paper I for the colour transformation of 
R and B literature values and the g measurements, and no 
correction between the R and r, and I and i band, respec- 
tively. Due to the fact that the photometric errors, of the 
input catalogues from Lauberts & Valentijn (198£) as well 



as within the RC3 catalogue (de Vaucouleurs et al. 1991) 
are of the order of 0.1 mag, we do not apply any further 
corrections. Galaxies calibrated in this way are marked 
with I in column (5) of Table ||, whereas for galaxies which 
did not have published values for their magnitude in the 
observed filter, we interpolated from calibrated images of 
the same night, by comparing the count rates for the sky 
value. These galaxies are marked with an i. For a few 
nights no galaxy with published photometry was observed 
and in these cases we used interpolated night sky values 
from the same observing run. Images calibrated in this way 
are marked with e in Table ^. The resulting zero points 
and central surface brightness values in these cases should 
be interpreted carefully, although the derived structural 
parameters, like scalelength and scaleheight, are not in- 
fluenced by any uncertainty in the flux calibration. 
Appendix A shows the contour plots and selected radial 
profiles for the 14 objects (25 images) not already pub- 
lished in Paper I. 



NASA/IPAC Extragalactic Database (NED) 
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Table 1. Global parameters of sample galaxies with models 



No. 


name 


filter 


tint. 


run 


RA 


DEC 


RC3 


T 


D25 


VQ 


D 








[min] 


ID 


(2000.0) 


type 




['] 


km s ^ 


Mpc 


(V 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


ESO 112-004 


r 


40 


E3* 


002804.2 


-580611 


.S.R6*. 


5.6 


1.32 


— 


— 


2 


ESO 150-014 


r 


20 


E3* 


003637.9 


-565424 


.L..+V 


-0.7 


1.91 


8257 


107.05 


3 


NGC 585 


R 


20 


LI 


013142.5 


-005555 


.S..1*/ 


1.0 


2.14 


5430 


72.15 


4 


ESO 244-048 


r 


15 


E3 


013908.8 


-470742 


.S..3./ 


3.0 


1.38 


6745 


87.09 


5 


NGC 973 


R 


10 


LI 


023420.2 


+323020 


.S..3.. 


3.0 


3.72 


4853 


66.28 


6 


UGC 3326 


R 


30 


LI 


053936.0 


+771800 


.S..6*. 


6.0 


3.55 


4085 


57.82 


7 


UGC 3425 


R 


30 


LI 


061442.0 


+663400 


.S..3.. 


3.0 


2.51 


4057 


57.04 


8 


NGC 2424 


R 


15 


LI 


074039.8 


+391359 


.SBR3*/ 


3.0 


3.80 


3113 


43.10 


9 


IC 2207 


R 


10 


L2 


074950.8 


+335743 


.S..6*. 


6.0 


2.04 


4793 


65.10 


10 


ESO 564-027 


r 


30 


E2 


091154.4 


-200703 


.S..6*/ 


6.3 


4.07 


2177 


26.93 


11 


ESO 436-034 


g 


60 


E3* 


103244.2 


-283646 


.S..3./ 


3.0 


2.09 


3624 


46.02 


12 


ESO 319-026 


g 


30 


E3* 


113020.0 


-410357 


.S..5./ 


5.3 


1.48 


3601 


45.32 


12 


ESO 319-026 


i 


30 


E3 
















12 


ESO 319-026 


r 


30 


E2 
















12 


ESO 319-026 


r 


30 


E3* 
















13 


ESO 321-010 


g 


30 


E3* 


121142.2 


-383253 


.S..0*/ 


0.0 


1.86 


3147 


39.53 


13 


ESO 321-010 


r 


30 


E3* 
















14 


NGC 4835A 


r 


40 


E3* 


125713.6 


-462243 


.S..6*/ 


6.0 


2.57 


3389 


42.54 


15 


ESO 575-059 


r 


15 


E2 


130744.5 


-192348 


.LA.+?/ 


-0.8 


1.86 


4570 


59.78 


16 


ESO 578-025 


g 


30 


E2 


140815.5 


-200019 


— 


1.0 


1.60 


6364 


83.93 


16 


ESO 578-025 


g 


30 


E3* 
















16 


ESO 578-025 


i 


30 


E2 
















16 


ESO 578-025 


r 


30 


E2 
















17 


ESO 446-018 


r 


30 


E2 


140838.7 


-293412 


.S..3./ 


3.0 


2.34 


4774 


62.15 


18 


IC 4393 


r 


30 


E2 


141749.5 


-312056 


.S..6?/ 


6.0 


2.40 


2754 


35.15 


19 


ESO 581-006 


r 


30 


E3* 


145803.1 


-192329 


.SBS7P/ 


7.0 


1.70 


3119 


40.91 


20 


ESO 583-008 


r 


30 


E3 


155750.5 


-222947 


.S? 


6.0 


1.51 


7399 


97.98 


21 


UGC 10535 


r 


25 


E2 


164600.0 


+062800 


.S..2.. 


2.0 


1.10 


7586 


102.48 


22 


NGC 6722 


r 


10 


E3* 


190339.6 


-645341 


.S..3./ 


3.0 


2.88 


5749 


73.79 


23 


ESO 461-006 


r 


60 


E3* 


195155.9 


-315852 


.S..5./ 


5.0 


1.62 


5949 


78.26 


24 


IC 4937 


g 


20 


El 


200518.0 


-561520 


.S..3./ 


3.0 


1.95 


2337 


28.68 


24 


IC 4937 


i 


20 


El 
















24 


IC 4937 


r 


30 


E3* 
















25 


ESO 528-017 


g 


30 


E3* 


203320.8 


-270549 


.SB.6?/ 


5.7 


1.59 


6115 


80.73 


25 


ESO 528-017 


i 


30 


E3 
















25 


ESO 528-017 


r 


60 


E3* 
















26 


ESO 187-008 


r 


30 


E3* 


204325.2 


-561217 


.S..6./ 


6.0 


1.51 


4412 


56.31 


27 


ESO 466-001 


i 


40 


E3 


214232.3 


-292210 


.S..2./ 


2.0 


1.38 


7068 


93.13 


28 


ESO 189-012 


g 


60 


E3* 


215538.7 


-545233 


.SA.5*/ 


5.0 


1.66 


8398 


109.04 


28 


ESO 189-012 


i 


20 


E3 
















28 


ESO 189-012 


r 


30 


E3* 
















29 


ESO 533-004 


r 


20 


El 


221403.2 


-265618 


.S..5*/ 


4.8 


2.34 


2594 


33.54 


30 


IC 5199 


g 


30 


E3* 


221933.0 


-373201 


.SA.3*/ 


3.0 


1.55 


5061 


65.78 


30 


IC 5199 


i 


30 


E3 
















31 


ESO 604-006 


r 


30 


E3* 


231454.0 


-205944 


.S..4./ 


4.0 


1.86 


7636 


100.92 



* marks images already published in Paper 1 



2.5. Distance estimates 

In order to derive the intrinsic values of the scale param- 
eters and to compare physical dimensions we tried to es- 
timate distances for our galaxies. Therefore we took pub- 
lished radial velocities corrected for the Virgo centric in- 



fall from the LEDA0 database, and estimated the distance 
following the Hubble relation with a Hubble constant of 
Hq = 75 km s-^Mpc-i. Table gives the heliocentric ra- 
dial velocities (11) according to LEDA, and our estimated 
distances (12). 

^ Lyon/Meudon Extragalactic Database (LEDA) 
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3. Disk models 

3.1. Background 

Our disk model is based upon the fundamental work of van 



der Kruit and Searle ( |l98l4 |l981b| , |19824 |l982b| ). They 
tried to find a fitting function for the light distribution in 
disks of edge-on galaxies. These galaxies are, compared to 
the face-on view, preferred for studying galactic disks due 
to the fact, that in this geometry it is possible to disentan- 
gle the radial and vertical stellar distribution. Their model 
include an exponential radial light distribution found for 
face-on galaxies (de Vaucouleurs, 1959; Freeman, 1970D , 
a sech^ behaviour in z, which is expected for an isother- 



mal po pulati on in a plan-parallel system (Camm, 195C ; 
Spitzer, 1942 ), an d a sh arp edge of the disk, first observed 
by van der Kruit ( 197E ) in radial profiles of edge-on galax- 
ies. The resulting luminosity density distribution for this 
symmetric disk model is (van der Kruit & Searle 1981a): 



L(R,z) — Lq exp ( — ^ I sech^f — 



R < Rco (1) 



L being the luminosity density in units of [Lq pc^^], Lq 
the central luminosity density, R and z are the radial resp. 
vertical axes in cylinder coordinates, h is the radial scale- 
length and Zq the scaleheight, and Rco is the cut-off radius. 
The empirically found exponential radial light distribution 
is now well accepted and it is proposed that viscous dissi- 
pation could be responsible (Firmani et al. 1996, Struck- 
Marcel |199lt Saio & Yoshu \199(\ , Lin & Pringle |1987D , 
although there is so far no unique explanation for the disk 
being exponential. An alternative description of the form 
1/R proposed by Seiden et al. ( |1984D did not get much 
attention, although it emphasizes the empirical nature of 
the exponential fitting function. 

To avoid the strong dust lane and to follow the light 
distribution down to t he reg ion z w Wainscoat (1986) 
and Wainscoat et al. (1989) carried out NIR observa- 



tions using the much lower extinction in this wavelength 
regime compared to the optical. They found a clear ex- 
cess over the isothermal distribution and proposed the z- 
distribution to be better fitted by an exponential func- 
tion f{z) = exp(— z/zo). According to van der Kruit 
( |1988| ) such a distribution would led to a sharp mini- 
mum of the velocity dispersion in the plane, which is 
not observed (Fuchs & Wielen 1987 ). Therefore he pro- 
posed /(z) — sech(z/zo) as an intermediate solution. De 
Grijs ( 1997| ) extended this to a family of density laws 
gmiz,zo) = 2-2/™ ffo se ch^/"' (77iz/2zo) (m > 0) follow- 
ing van der Kruit ( |1988| ), where the isothermal (m = 1), 
and the exponential (m = oo) cases represent the two ex- 
tremes. 

Therefore the luminosity density distribution can be writ- 
ten as: 



L{R, z) — Lq exp 



with H(a;o — x) being the Heaviside function. 
In order to limit the choice of parameters we restrict 
our models to the three main density laws for the z- 
distribution (exponential, sech, and sech^). Due to the 
choice of our normalised isothermal case zq is equal to 2hz , 
where is the usual exponential vertical scale height: 

= 4 exp (-2 ^ 
V ^0 

, / 2z 
~ 2 sech I — 
Zq 



fiiz 

h{z) = sech^l — 



In contrast to Paper I and Barteldrees & Dettmar ( |1989|) 
we define the cut-off radius at the position where the ra- 
dial profiles become nearly vertical, corresponding to the 
mathematical description. They tried to avoid any confu- 
sion due to the lower signal-to-noise in the outer parts, by 
fixing the cut-off radius where the measured radial profile 
begin to deviate significantly from the pure exponential 
fit. 

3.2. Numerical realisation 

The model of the two dimensional surface photometric in- 
tensity results from an integration along the line of sight of 
the three dimensional luminosity density distribution ^ 
with regard to the inclination i of the galaxy. Describing 
the luminosity density of the disk in a cartesian-coordinate 
grid K(x-t/-z) with R = {x^ + y^) leads to the follow- 
ing coordinate transformation into the observed inclined 
system K'(a;'-y'-z') with x' pointing towards the observer, 
whereas the rotation angle between the two systems cor- 
responds to 90° — i: 

X — x' sin(i) — z' cos(i) 
2/ = y' 

z = cos(i) -I- z' sin(i) 

Taking into account this transformation we have to inte- 
grate equation (H), obtaining an equation for the intensity 
of the model disk depending on the observed radial and 
vertical axes y' and z' on the CCD: 



+ 00 



l{v\z')^ / L{x{x',z',i),y',z{x',z',i))dx' 



(3) 



Together with equation (^ this gives: 



y/{x' sini-z' cos?)2+y'2\ 
Lq I exp I /„(z , zq) 



/„(z,zo) R{Rco~R) 



(2) 



■■ H ( i?co - \/ix' sini - z' cosi)^ + y'^ ) dx' (4) 
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Therefore six free parameters fit the observed surface in- 
tensity on the chip (y', z' plane) to the model: 

/ = /(y', z', I, n, Xo, -Rco, zq) (5) 

Figure shows a sequence of computed models with an 
isothermal z-distribution (n = 3) for the exact edge-on case 
[i — 90°) with characteristic values for the ratio Rco/h: 
1.40,2.88,5.00 and for h/zo: 2.0,4.0,7.3 keeping the cut- 
off radius i?co and the total luminosity Ltot oc zh^Lo con- 
stant. The latter is causing a different central surface 
brightness fiQ for each model, ranging from 23.8 to 19.3 
starting with /io ~ 21.2 mag/D" for the reference model 
with Rco/h = 2.88 and h/zo = A.O. All contour lines falling 
within the interval fi = 25.0 - 19.5 are plotted with a spac- 
ing of 0.5. 

3.3. Method 1 

3.3.1. Determination of fitting area 

The first step is to divide the galaxy into its quadrants. 
The four images are then averaged according to their 
orientation, follo wing van der Kruit & Searle ( |l981aD 
and Shaw ( 1993| ). Thereby larger foreground stars and 
asymmetrical perturbations in the intensity distribution 
are eliminated by omitting this region during averaging. 
Smaller foreground stars are removed by median filtering. 
The average quadrant should result at least from three 
quadrants to get a representative image of the galaxy. This 
averaging additionally increases the signal-to-noise ratio. 
In order to determine the fitting area for modelling the 
disk component on the final quadrant, one has to avoid 
the disturbing influence of the bulge component and the 
dust lane. The region dominated by the bulge is fixed fol- 
lowing Wyse et al. (1997) defining the bulge component by 
"light that is in excess of an inward extrapolation of a con- 
stant scale-length exponential disk" . Therefore the clear 
increase of the intensity towards the center which can be 
seen in radial cuts determines an inner fitting boundary 
-Rmin - We tried to minimize the dust influence (cf. Sec- 
tion 4.4.2 ) by placing a lower limit z^in in the vertical 
direction by visual inspection. Additionally we restricted 
the remaining image by a limiting contour line /iiim, where 
the intensity drops below a limit of 3cr on the background. 
We are aware of the problem that these are rather rough 
definitions difficult to reproduce without quoting the ex- 
act values for i^min, ^min, Miim- Howevcr, the final choice of 
the fitting area is a complex and subjective procedure de- 
pending on the intrinsic shape of each individual galaxy, 
the influence of their environment, and the quality of the 
image itself. Therefore it is not possible to quote exact 
general selection criteria and to derive the structural pa- 
rameters straight forward. One solution is to do it in a 
consistent way for a large sample, leaving the problem of 
comparing results from different methods (cf . Section 



3.3.2. Numerical fltting 

The numerical realisation of the fltting procedure mini- 
mizes the difference {SQ) between the averaged quadrant 
and a modelled quadrant based on equation (W). 



SQ 



E 

3 



log {Io, {yj,Zj)) - log {Im, {yj,Zj)) 



loj is the intensity within the average quadrant (observed 
intensity) and Im of the modelled intensity. In contrast 
to Shaw & Gilmore ( |1989| , |1990D and Shaw ( |1993D using 
a similar approach for their models we do not weight in- 
dividual pixels. They weight the difference by the error 



in the surface brightness measure, which Shaw (1993) de- 
rives from the averaging of the quadrants. However, this 
method implies an absolute symmetry for the disk in z 
and y, which is not the case for real galaxies. These kind 
of errors only reflect the asymmetry of the galaxy. Using 
the observed errors in the surface brightness for weighting 
individual pixels does not result in a considerable advan- 
tage because they are nearly the same after smoothing. 
The minimal SQ is found by varying five of the six free 
parameters of the model (cf. eq. whereas the parameter 
Rco is determinated by cuts parallel to the major axis. A 
significant decrease of the intensity extrapolated t o / = 
gives the value of i?co (van der Kruit & Searle, 1981a ), 
therefore it is important that the intensity at the cut-off 
radius is well above the noise limit. The other five pa- 
rameters are determined by fixing the smallest SQ. For 
the three different functions f{z) and every possible i 
{Ai = 0.5°) the remaining parameters {Lq, h, and zq) 
are varied with the "downhill simplex-method" (Nelder & 
Mead, 1965| ; Press et al., 1988) until the global minimum 
of SQ is found. is calculated by a numerical gaussian 
integration of equation (^. The possible inclination angles 
can be restricted from the dust lane of the galaxy (Paper 
I). Tests of model disks with additional noise show that 
the "downhill simplex-method" found the input inclina- 
tion i, the used model /(z) and the other disk parameters 
within errors of (5/i, Szq, and SLq < 1%. An estimation of 
the errors of the parameters for the best model disk can be 
made by inspection of the parameter space for SQ around 
the smallest SQ, with slightly different fitting areas, and 
different values of i?co- f{z) is in almost all cases the same 
and the variation of i is only small (±1°). The differences 
in h and Zq are about 15% (in some cases up to 25%), and 
Lq varies about a factor of 2. 

3.4. Method 2 

Within this method all disk parameters as well as the 
choice of the optimum function /(z) are estimated by a 
direct comparison of calculated and observed disk profiles 
by eye. 

As a first step, the inclination i is determined by using the 
axis ratio of the dust lane. Depending on the shape of the 
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Fig. 1. Different galaxy models computed according to equation 
columns, or within the rows, whereas Ltot and i?co are constant for all models using an isothermal description for 
the exact edge-on case 
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major axis (average) at z=2.376 minor axis (average) at R=54.72 
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Fig. 2. Example of the fitting procedure: Set of typical 
radial (left panel) and vertical (right panel) profiles used 
for the fitting process of the galaxy ESO 461-006 together 
with the model. 



dust lane it is possible to restrict the inclination to ±1.5°. 
The central luminosity density Lq is calculated automat- 
ically for each new parameter set by using a number of 
preselected reference points along the disk. Given a suf- 
ficient signal-to-noise ratio, the cut-off radius Rco can be 
determined from the major axis profile. Thus, the remain- 
ing fitting parameters are the disk scale length and height, 
h and zq, as well as the set of 3 functions /(z). The scale 
length h is fitted to a number (usually between 4 and 6) of 
major axes profiles (left panel of Fig. ^). The quality of the 



fit for the vertical profiles along the disk can be used as a 
cross-check for the estimated scalelengths (right panel of 
Fig. H) . The disk scaleheight zq is estimated by fitting the 
z-profiles, outside a possible bulge or bar contamination. 
The vertical disk profiles of most of the galaxies inves- 
tigated enable a reliable choice of the quantitatively best 
fitting function f{z). This is due to the fact that the devia- 
tions between different functions become visible at vertical 
distances larger than that of the most sharply-peaked dust 
regions. The first raw fitting steps are usually carried out 
by using a reduced number of both major and minor axis 
profiles simultaneously. Afterwards, when a good first fit- 
ting quality is reached, a complete set of major and minor 
calculated and observed axes profiles are investigated in 
detail. At the end of this procedure a final, complete disk 
model is calculated using all previously estimated disk pa- 
rameters. 

4. Results 

4-.1. Distribution of disk parameters 

Table || contains the best fit model for each image. To- 
gether with the galaxy name (i), the filter (2), and the 
referring image (3), with integration time, and run ID, we 
list the inclination (4), the best fitting function for the 
z-distribution (5), the calibration index (6) (ref. Section 
^.4[ ), and the central surface brightness of the model (7), 
without correcting for inclination. According to the dis- 
tance tabulated in Table |^, the cut-off radius -Rco (s) is 
given in kpc and arcsec as well as the scalelength h (9) 
and the vertical scaleheight zq (10) which is normalised 
to the isothermal case, being two times an exponential 
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Table 2. Determined parameters set 



galaxy 


filter 


image 


i 




cali. 


Ho 


Rco 




h 




zo 








n 






mag/n" 


["] 


[kpc] 


["] 


[kpc] 


["] 


[kpc] 


(1) 


(2) 


(3) 


U) 


(5) 


(6) 


CJ) 


(S) 




(9) 


(10) 


ESO 112-004 


r 


40E3 


87.5 


sech 


1 


21.14 


45.0 


— 


22.7 


— 


2.8 


— 


ESO 150-014 


r 


20E3 


90.0 


sech 


1 


22.00 


64.1 


33.27 


23.4 


12.14 


5.5 


2.85 


NGC 585 


R 


20L1 


88.0 


sech 


e 


21.76 


70.0 


24.49 


32.7 


11.44 


10.6 


3.71 


ESO 244-048 


r 


15E3 


87.0 


exp 


1 


20.54 


45.4 


19.17 


13.7 


5.78 


7.2 


3.04 


NGC 973 


R 


lOLl 


89.5 


exp 


e 


20.83 


105.0 


33.74 


51.5 


16.55 


12.4 


3.99 


UGC 3326 


R 


30L1 


88.0 


sech^ 


e 


21.50 


101.5 


28.45 


70.4 


19.74 


4.7 


1.32 


UGC 3425 


R 


30L1 


87.0 


sech^ 


e 


21.01 


80.5 


22.26 


29.2 


8.08 


7.4 


2.05 


NGC 2424 


R 


15L1 


86.5 


exp 


e 


20.52 


112.0 


23.40 


31.7 


6.62 


11.9 


2.49 


IC 2207 


R 


10L2 


86.5 


exp 


e 


21.04 


56.0 


17.67 


38.2 


12.06 


6.7 


2.12 


ESO 564-027 


r 


30E2 


88.0 


sech 


1 


21.11 


140.4 


18.33 


50.9 


6.65 


6.6 


0.86 


ESO 436-034 


g 


60E3 


88.0 


sech 


1 


20.96 


82.1 


18.32 


22.6 


5.04 


6.8 


1.52 


ESO 319-026 


g 


30E3 


86.5 


sech 


e 


21.51 


60.1 


13.21 


14.2 


3.12 


3.1 


0.68 


ESO 319-026 


r 


30E2 


86.0 


exp 


i 


21.14 


63.0 


13.84 


14.7 


3.23 


3.2 


0.70 


ESO 319-026 


r 


30E3 


88.0 


sech^ 


i 


21.94 


61.9 


13.60 


13.2 


2.90 


2.8 


0.62 


ESO 319-026 


i 


30E3 


88.0 


sech 


1 


21.01 


64.8 


14.24 


14.3 


3.14 


3.4 


0.75 


ESO 321-010 


g 


30E3 


88.0 


sech 


1 


20.11 


64A 


12.34 


19.9 


3.81 


4.7 


0.90 


ESO 321-010 


r 


30E3 


88.0 


sech 


1 


19.54 


64.8 


12.42 


21.6 


4.14 


5.0 


0.96 


NGC 4835A 


r 


40E3 


85.5 


exp 


1 


20.92 


90.0 


18.56 


40.4 


8.33 


7.9 


1.63 


ESO 575-059 


r 


15E2 


87.0 


sech^ 


1 


20.78 


60.5 


17.53 


22.6 


6.55 


6.4 


1.86 


ESO 578-025 


g 


30E2 


86.5 


exp 


1 


21.58 


50.7 


20.63 


17.0 


6.92 


7.5 


3.05 


ESO 578-025 


g 


30E3 


86.5 


exp 


1 


21.69 


50.7 


20.63 


17.0 


6.92 


7.5 


3.05 


ESO 578-025 


r 


30E2 


86.0 


exp 


1 


21.02 


50.4 


20.51 


14.9 


6.06 


7.4 


3.01 


ESO 578-025 


i 


30E2 


86.0 


exp 


1 


20.27 


47.5 


19.33 


20.3 


8.26 


7.6 


3.09 


ESO 446-018 


r 


30E2 


86.5 


sech 


1 


20.45 


75.6 


22.78 


23.7 


7.14 


3.9 


1.18 


IC 4393 


r 


30E2 


87.0 


exp 


1 


20.30 


75.6 


12.88 


29.8 


5.08 


5.8 


0.99 


ESO 581-006 


r 


30E3 


86.5 


exp 


1 


21.33 


55.8 


11.07 


18.2 


3.61 


5.6 


1.11 


ESO 583-008 


r 


30E3 


87.0 


exp 




20.84 


56.2 


26.70 


14.0 


6.65 


3.8 


1.81 


UGC 10535 


r 


25E2 


88.0 


sech'^ 




21.50 


41.4 


20.57 


11.4 


5.66 


4.3 


2.14 


NGC 6722 


r 


10E3 


86.5 


sech 


1 


19.47 


86.0 


30.77 


21.6 


7.73 


6.2 


2.22 


ESO 461-006 


r 


60E3 


87.5 


exp 


1 


20.96 


61.6 


23.37 


21.1 


8.01 


3.8 


1.44 


IC 4937 


g 


20E1 


88.5 


soch'^ 


1 


22.73 


74.9 


10.41 


93.8 


13.04 


5.6 


0.78 


IC 4937 


r 


30E3 


88.0 


sech^ 


1 


21.29 


83.5 


11.61 


32.1 


4.46 


5.9 


0.82 


IC 4937 


i 


20E1 


88.5 


sech 




19.92 


78.1 


10.86 


38.6 


5.37 


7.2 


1.00 


ESO 528-017 


g 


30E3 


86.5 


exp 


1 


21.98 


57.6 


22.54 


20.5 


8.02 


3.5 


1.37 


ESO 528-017 


r 


60E3 


86.5 


exp 


1 


21.28 


55.1 


21.56 


19.9 


7.79 


2.7 


1.06 


ESO 528-017 


i 


30E3 


86.5 


sech 




20.89 


50.0 


19.57 


22.6 


8.85 


3.2 


1.25 


ESO 187-008 


r 


30E3 


85.5 


exp 




20.77 


50.0 


13.65 


15.1 


4.12 


4.4 


1.20 


ESO 466-001 


i 


40E3 


87.0 


exp 




19.52 


52.6 


23.76 


13.0 


5.87 


8.2 


3.70 


ESO 189-012 


g 


60E3 


87.0 


exp 




21.50 


56.2 


29.81 


26.8 


14.21 


3.8 


2.02 


ESO 189-012 


r 


30E3 


86.5 


exp 




20.71 


56.5 


29.97 


20.6 


10.93 


3.4 


1.80 


ESO 189-012 


i 


20E3 


87.0 


sech 




20.44 


54.7 


29.01 


22.7 


12.04 


3.2 


1.70 


ESO 533-004 


r 


20E1 


88.0 


exp 




20.32 


68.4 


11.12 


33.1 


5.38 


7.7 


1.25 


IC 5199 


g 


30E3 


86.5 


exp 




21.55 


64.1 


20.44 


19.9 


6.35 


5.6 


1.79 


IC 5199 


i 


30E3 


86.5 


exp 




19.97 


57.6 


18.37 


19.9 


6.35 


4.9 


1.56 


ESO 604-006 


r 


30E3 


90.0 


sech 




21.22 


70.6 


34.54 


27.9 


13.65 


3.8 


1.86 



scaleheight hz- For the seven galaxies with available im- 
ages in more than one filter, wc do not sec any correlation 
of fitted parameters with different wavelength, although 
we find the same inclination angle for the best fitted disk 
within the range of the errors. Appendix A shows the best 
fitting model as an overlay to selected radial profiles for 
each image. The subsequent analysis of the distribution 



for the different parameters concerning the formation and 
evolution of galaxies will be given in forthcoming papers. 

4-2. Comparison of different methods 

The different fitting methods were independently devel- 
oped within two diploma theses (Liitticke 1996, Schwarz- 
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kopf 1996). The quality of the data basis for each project 
was the same. From the sample presented here there were 
five objects in common. These are used to compare the 
two methods and determine the quantitative difference of 
the derived parameters. Table || shows the results for the 

Table 3. Comparison of the different determined param- 
eter sets for the same galaxy images 



definition of the cut-off radius. 

Only three galaxies are in common with the sample of de 



galaxy f{z) 


i 


z 


h 


Rco 






["] 


["1 


["] 



ESO 321-010 



ESO 461-006 r 



IC 4937 r 



ESO 189-012 



ESO 604-006 



sech 
exp 
exp 
exp 
exp 
sech 
exp 
exp 
sech 
sech'^ 



88.0 
88.0 
87.5 
87.5 
86.5 
89.5 
86.5 
88.0 
88.5 
89.5 



5.0 
5.2 
3.8 
3.6 
6.1 
7.0 
3.2 
3.6 
4.1 
3.2 



21.6 
25.9 
21.1 
16.2 
26.8 
27.4 
19.3 
21.3 
39.8 
24.5 



64.8 
64.1 
61.6 
61.2 
83.5 
75.6 
56.2 
56.9 
67.7 
73.4 



five images. The mean deviation in the determined incli- 
nation is « 1° and 12.4% for the scaleheight (ranging from 
5.0%-26.6%) whereas for three images different functions 
for the z distribution were used. The mean difference for 
the radial scalelength is 20.6% (2.1%-47.2%) and 4.2% for 
the determination of the cut-off radius. 
A subsequent analysis shows that it is not possible to as- 
cribe the sometimes quite large discrepancies to the qual- 
ity of the individual method. It turns out, that the main 
problem is the non-uniform determination of the fitting 
area. The intrinsic asymmetric variations of a real galactic 
disk compared to the model enforce a more subjective re- 
striction of the galaxy image to the fitting region, whereby 
for example one has to exclude the bulge area and the dust 
lane. 

This finding is in agr eement with the study of Knapen & 
van der Kruit (1991) who compared published values of 
the scalelength and find an average value of 23% for the 
discrepancy between different sources. As already men- 
tioned by Schombert & Bothun ( 1987 ) the limiting factor 
for accuracy of the decomposition is not the typical S/N 
from the CCD-telescope combination nor the errors in the 
determination of the sky background, but the deviation of 
real galaxies from the standard model. 



4-.3. Comparison with the literature 

In our former study (Paper I) with an earlier method to 
adapt equation (|^), 20 of our 45 galaxy images have al- 
ready been used. We decided to re-use them in this study 
to get models for as many galaxies as possible in a homo- 
geneous way. Additionally, Paper I only presents the best 
fit values for the isothermal model, and uses a different 



Grijs ( |1998| ): ESO 564-027, ESO 321-010, and ESO 446- 
018. The mean difference for the scalelength is 10.4 % 
(ranging from 1.2%-20.3%) and for the scale height (nor- 
mafised to the isothermal case) 4.0% (0.0%-8.2%). For the 
remaining galaxies there are no models in the literature. 

4- 4- Model limitations 

Our model only represents a rather simple axisymmetric 
three-dimensional model for a galactic disk, consisting of 
an one component radial exponential disk with three dif- 
ferent laws for the density distribution in the z-direction 
and a sharp outer truncation. Therefore it does not include 
additional components, such as bulges, bars, thick disks, 
or rings, and cannot deal with any asymmetries. Features 
like spiral structure or warps are not included, whereas 
Reshetnikov & Combes (1998) multiply their exponential 
disk by a spiral function introducing an expression to char- 
acterise an intrinsic warp depending on the position angle 
outside of a critical radius. 

The choice of our fitting area tries to avoid the dust lane, 
possible only for almost edge-on galaxies, as a firs t step 
to account for the dust influence (cf. Section 4.4.2| ). Ex- 
amples of models including a radiative transfer with an 
extinction coefficient kx{R, z) can be found in Xilouris et 
al. (1999). However, introducing more and more new com- 
ponents and features automatically increases the amount 
of free parameters. Therefore we restricted our model to 
the described six parameters, to obtain statistically mean- 
ingful characteristics for galactic disks. 
In the following we demonstrate that a simple disk model 
omitting the bulge component and the dust lane give in- 
deed reasonable parameters. 



4.4.1. The influence of the bulge component 

We have studied the influence of the bulge for some of 
our objects including the earliest type galaxy in our sam- 
ple (ESO 575-059) presented here. We have subtracted 
our derived disk model from the galaxy and then tried to 
find the best representation for the remaining bulge by a 
de Vaucouleurs r^/^ or an exponential model. Taking the 
slope of the vertical profile at i? = and a fixed axis ra- 
tio, we have constructed the 2-dimensional model of the 



bulge. In agreement with Andredakis et al. (1995) we find, 
that bulges of early type galaxies are better fitted by an 
exponential profile than by a r^^^. Figure ^ shows the re- 
sulting vertical and radial cuts for ESO 575-059 together 
with the models. Despite the deviation between R = 10" 
and R = 20" which could be attributed to an additional 
component (inner disk or bar), we do not find any evi- 
dence for changing our disk model due to the influence of 
the bulge. 

Therefore we conclude, that it is possible to nearly avoid 
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Fig. 3. Minor and major axes profile for ESO575059rl5E2 
(solid lines) together with the best disk model (dotted), 
the best bulge model (dashed), and the resulting combined 
profile (dotted-dashed). 



any influence of the central component by fitting outside 
the clearly visible bulge region. 

4.4.2. The influence of dust 

Dust disturbs the light profile by a combination of absorp- 
tion and extinction and the net effect has to be calculated 
by radiation transfer models. Therefore it is not obvious 
that outside the "visible" dust lane, which is excluded for 
the fitting area as a first step, the dust will not play a 
major role in shaping the light distribution. Xilouris et al. 
( |1999D Bian chi et al. ( |1996| ), de Jong ( |l996aD , and Byun 
et al. ( |1994 ) have recently addressed this problem in more 
detail. Although they investigated the influence of the dust 
on the light distribution by quoting best fit structural pa- 
rameters for the star-disk as well as the dust-disk, they 
did not quantify the influence on the star-disk parameters 
derived by standard fitting methods without dust. Even 
Kylafis 



Bahcall (1987) state within their fundamental 
paper on finding the dust distribution for NGC 891 that 
"in order to avoid duplication of previous work, we will 
take ... (the values for the star-distribution estimated by 



standard fitting methods)" . 

We checked the influence of the dust on our determined 
parameter set, by studying simulated galaxy images with 
three different dust distributions. These dusty galaxies 
were kindly provided by Simone Bianchi who calculated 
images with known input parameters for the star and 
dust distribution with his Monte Carlo radiative transfer 



method (Bianchi et al. 1996). We have defined a worst, 



of our mean 
The worst 



best, and transparent case, according to t he dust distri- 
butions presented by Xilouris et al. ( 1999| ) 
stellar disk (i?co = 2.9, h/z = 4, and fi{z)) 
case is calculated with = 0.51, h^/h^, — 1.55, and 
Zd/ z^, = 0.75, the best case with tr = 0.20, hd/h^, — 1.08, 
and ZdZ-z* = 0.32, and a transparent case without dust. 
To be comparable we used the same method for selecting 
the fitting region by masking the "visible" dust lane and 
reserving a typical area for a possible bulge component 
using mean values for the transparent case. In contrast to 
our standard procedure we do not restrict the inclination 
range from the appearance of the dust lane, but specify 
the best SQ model in the range i = 80°-90° in Table |. For 
the models marked with a we pretend the correct input 
inclination. Table ^ demonstrates, that even for the worst 
case we are able to reproduce the input parameters within 



the range of the typical 20% error discussed in Section 4.2 
It should be mentioned that for each case we overestimate 
the input scalelength and -height, whereas the determina- 
tion of Rco does not depend on the dust distribution. The 
implication on the distribution of the ratio Rco/h will be 
discussed in a forthcoming paper. 

Table 4. Comparison of the results for the stellar disk 
(input: =8.4, /i, =33.7, and /(z)=exp) from our fitting 
procedure for the three different dust distributions. 



case 


iin 


fn 


^out 




ft. 




< 




[°] 




[°] 


["] 


["] 


[%] 


[%] 


trans. 


87.5 


1 


87.5 


8.5 


34.0 


1.2 


0.9 


best 


85.0 


1 


84.5 


8.5 


35.0 


1.2 


3.9 


best 


87.5 


1 


86.5 


8.6 


35.5 


2.4 


5.3 


best 


90.0 


1 


87.0 


8.7 


36.0 


3.6 


6.8 


best-k 


90.0 


1 


90.0 


8.8 


35.0 


4.8 


3.9 


worst 


85.0 


2 


85.0 


8.9 


39.9 


6.0 


18.4 


worst 


87.5 


2 


86.0 


8.8 


39.1 


4.8 


16.0 


worst 


90.0 


2 


88.0 


9.7 


39.4 


15.5 


16.9 


worst* 


90.0 


2 


90.0 


9.8 


39.1 


16.7 


16.0 



4-5. Comments on individual galaxies 

Trying to adapt a simple, perfect, and exact symmetric 
model to real galaxies always implies a compromise be- 
tween the degree of any deviation and the final model 



(Section 4.6). The following list will provide some typical 
caveats found during the fit procedure which will char- 
acterise the quality of the specified model for individual 
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galaxies. 

ESO 112-004: warped, asymmetric, central part slightly 
tilted compared to disk, after fitting still remaining resid- 
uals 

ESO 150-014: slightly warped, minor flatfield problems 
NGC 585: remaining residuals 

ESO 244-048: possible two component system, slope of 
inner radial profile significantly higher than of an outer 
one, final model fits the inner parts 

NGC 973: one side disturbed by stray light of nearby 
star, seems to be radially asymmetric, remaining residu- 
als 

UGC 3425: superimposed star on one edge 

NGC 2424: model does not fit very well without obvious 

reason 

ESO 436-034: strong bulge component, possibly barred, 
hard to pinpoint final model, remaining residuals 
ESO 319-026: outer parts show u-shaped behaviour, re- 
maining residuals, therefore large (±2°) difference in in- 
clination angle 

ESO 321-010: u-shaped, no clear major axis visible, 
therefore uncertain rotation angle, bar visible, bulge ro- 
tated against disk 
NGC 4835A: strong residuals 

ESO 446-018: the different sides of the disk are asym- 
metric visible in radial profiles and on the contour plot 
IC 4393: similar to NGC 4835A 

ESO 581-006: galaxy shows typical late type profile, R^o 
questionable, but nevertheless final model seems to fit well 
ESO 583-008: disturbed by superimposed star, shows 
warp feature and a bar structure, i?co questionable, re- 
maining residuals 

UGC 10535: one side slightly extended 
NGC 6722: only one side observed, bulge rotated against 
disk, barred, strongly disturbed by dust absorption, ra- 
dial extension visible, therefore Rco should be treated with 
caution 

ESO 461-006: minor flatfield problem seems to cause 
asymmetry, although model looks fine 
IC 4937: similar to NGC 6722, dominating bulge, small 
disk, model significantly different compared to the i and r 
image, model possibly hampered by strong dust lane 
ESO 578-025: bar visible 

ESO 466-001: maybe two components, final model rep- 
resents only inner part, outer part clearly different from 
normal disk component 
ESO 189-012: slightly warped 

ESO 533-004: similar to NGC 4835A, model fits the 
whole galaxy, leaving more or less no bulge component 
IC 5199: slightly radial asymmetric 
ESO 604-006: only one side observed, bar structure vis- 
ible 



4-6. Comments on some rejected galaxies 

The model limitations described above constrain the ap- 
plication of our fitting process. Therefore we had to ex- 
clude about 20 galaxies from our original sample. They 
all show significant deviations from the simple geometry 
and an inclusion of their parameters obtained by forcing 
the model to fit the data will spoil the resulting parameter 
distribution. 

One larger group classified mainly as SO galaxies (e.g. 
NGC 2549, ESO 376-009, NGC 7332, ESO 383-085, ESO 
506-033) shows a completely different behaviour of the lu- 
minosity distribution in the outer parts compared to the 
other galaxies. They all show an additional component, 
mainly characterised as an elliptical envelope. This is al- 
ready visible in the contour plot, but becomes even more 
evident in a radial cut parallel to the major axis. In these 
cases the usual common curved decline of the profile (e.g. 
ESO 578-025) is missing, and is replaced by a more or less 
straight decline into the noise level, sometimes even by an 
upwards curved profile. Fitting these luminosity distribu- 
tion by our one component exponential disk with cut-off, 
will therefore naturally provide parameters qualitatively 
different compared to late type disks. This will be dis- 
cussed in detail in a forthcoming paper. 
Another group consists of galaxies dominated mainly by 
their bulges, whereas the disk is only an underlying com- 
ponent, partly characterized as having thick boxy bulges 
(Dettmar & Liitticke |f999| ), e.g. IC 4745, ESO 383-005, 
although there are also pure eUiptical bulges (e.g. ESO 
445-049, NGC 6948). 

In the case of ESO 383-048 and ESO 510-074 the radial 
profiles clearly indicate that a more complex model will 
be needed to fit these kind of multicomponent galaxies. 
Galaxies hke UGC 7170 or ESO 113-006 were excluded 
due to their strong warps, which made it impossible to fit 
the model in a consistent way. Mainly late type galaxies 
such as ESO 385-008, IC487I, UGC I28I, or ESO 376-023 
show a too patchy and asymmetric light distribution, that 
any attempt to fit the profiles will give only very crude, 
low quality parameters. UGC 11859 and UGC 12423 were 
rejected due to their thin faint disks, which will maybe 
overcome by taking new images with longer integration 
times to get a higher signal-to-noise ratio, whereas NGC 
5 193 A is completely embedded into the surface brightness 
distribution of its near companion. 
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Appendix A: Contour plots and radial profiles 

The following figures show in the top panel selected cuts 
parallel to the major axis (full line) as well as the best 
fit model (dashed line) according to Table ||. The bottom 
panel shows isophote maps of the surface brightness for the 
galaxies, rotated to the major axis. The faintest plotted 
contour is defined by a 3(t criterion of the background. 



Table Al lists the name (i), image (2), and the vertical 
positions of the plotted profiles parallel to the major axes 
in arcsec (3). The contours are plotted equally spaced with 
0.5 mag/n" starting from the limiting surface brightness 

Mlim. U)- 



Table Al. Plotted parameters 



galaxy image z /inm 

[ " ] [mag/D"] 
(1) (Z) (3) (4) 



NGC 585 


R20L1 


ESQ 244-048 


rl5E3 


NGC 973 


RlOLl 


UGC 3326 


R30L1 


UGC 3425 


R30L1 


NGC 2424 


R15L1 


IC 2207 


R10L2 


ESQ 564-027 


r30E2 


ESQ 319-026 


i30E3 


ESQ 319-026 


r30E2 


ESQ 575-059 


rl5E2 


ESQ 578-025 


g30E2 


ESQ 578-025 


i30E2 


ESQ 578-025 


r30E2 


ESQ 446-018 


r30E2 


IC 4393 


r30E2 


ESQ 583-008 


r30E3 


UGC 10535 


r25E2 


IC 4937 


g20El 


IC 4937 


i20El 


ESQ 528-017 


i30E3 


ESQ 466-001 


i40E3 


ESQ 189-012 


i20E3 


ESQ 533-004 


r20El 


IC 5199 


i30E3 



3,6,9,12 


24.7 


-3,-9 


24.8 


-2,-4,-6,-8 


23.8 


3,6,9 


24.8 


-1,-5,-9 


24.2 


1,3,5,7 


24.5 


2,6,10 


24.1 


2,4,6 


25.2 


2,4,6 


24.4 


2,4,6 


24.7 


3,6,9 


24.4 


2,4,6 


25.2 


2,4,6 


23.3 


2,4,6 


24.5 


2,4,6,8 


25.1 


-3,±6,±9 


25.1 


2,4,6 


25.0 


2,4,6 


25.5 


-2,-4,-6 


25.1 


-2,-4,-6 


23.3 


2,4,6 


24.2 


1,3,5 


23.8 


-1,-3,-5 


24.1 


-3,-6,-9 


24.6 


±2,±4,±6 


23.8 



remaining images already published in Paper I 
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